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Quantum-enhanced measurements use quantum mechanical effects in order to enhance the sensi- 
tivity of the measurement of classical quantities, such as the length of an optical cavity. The major 
goal is to beat the standard quantum limit (SQL), i.e. an uncertainty of order 1/^/N, where N is 
the number of quantum resources (e.g. the number of photons or atoms used), and to achieve a 
scaling 1/A^, known as the Heisenberg limit. So far very few experiments have demonstrated an 
improvement over the SQL. The required quantum states are generally highly entangled, difficult 
to produce, and very prone to decoherence. Here, we show that Heisenberg-limited measurements 
can be achieved without the use of entangled states by coupling the quantum resources to a com- 
j^ mon environment that can be measured at least in part. The method is robust under decoherence, 

and in fact the parameter dependence of collective decoherence itself can be used to reach a 1/iV scaling. 
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^\ ] Quantum mechanical noise imposes fundamental limitations on any measurement. The best-known example is the 

Heisenberg uncertainty relation, which provides a lower bound on the product of fluctuations of two non-commuting 

quantized variables. But even a classical system parameter x can in general not be measured with arbitrary precision 

p^ ; with a finite number of measurements due to the statistical nature of any quantum state. A lower bound on the 

I ' uncertainty is given by the smallest Sx such that two quantum states p{x) and p{x-\-Sx) lead to statistically significant 

(— I . differences for an optimally chosen observable. A similar problem exists in classical statistical analysis, where one 

d ' wants to distinguish between two probability distributions P{x) and P{x + dx), and the celebrated Cramer- Rao bound 

^ ] sets an ultimate lower bound on the uncertainty of a measurement of x based on the distinguishability of P{x) and 

, '^, '' P{x + Sx) [1[. That analysis has been generalized to the quantum world and has become known as "quantum 

[ parameter estimation theory" . Recently, the theory was used to prove lower bounds on the smallest measurable 

^s) . Sx for unitary time evolution. It was shown that if N replicas of the quantum system evolve independently and 

^ ' linearly, for an initially separable state no uncertainty smaller than l/yN can be achieved, i.e. the SQL, no matter 

C^ , how sophisticated the measurement. If initially entangled states are allowed, a 1/N scaling of the uncertainty is 

"^ ' the ultimate lower bound under otherwise identical conditions 0-ly]. The use of non-classical states of light for 

^f I Heisenberg hmited interferometry, notably the use of squeezed states, was proposed theoretically already in 1981 Q. 

So-called NOON states have been investigated for super-resolution |8l-[ig|. However, decoherence of these highly non- 



*^ I classical states has so far prevented reaching an uncertainty that scales as 1/N for systems with TV ^ 1 [Tl|, \l^- In 
pZ. • |13| entanglement-free Heisenberg-limited sensitivity of a phase-shift measurement was reported for several hundred 
^"'^ ' quantum resources by passing light many times through the phase-shifter. 

Decoherence arises when a quantum system interacts with an environment with many uncontrolled degrees of 

freedom, such as the modes of the electromagnetic field, phonons in a solid, or simply a measurement instrument |14| . 

k> , Decoherence destroys quantum mechanical coherence, and plays an important role in the transition from quantum 

5_j ■ to classical mechanics [15|. It becomes extremely fast for a mesoscopic or even macroscopic "distance" between the 

Cu ] components of a "Schrodinger cat" -t ype sup erposition of quantum states. Universal power laws rule the scaling of 

the decoherence rates in this regime |16l . 1 1711. Only recently could the collapse be time-resolved in experiments with 

relatively small "Schrodinger cat" -states |18l. Il9l|. However, decoherence can depend very sensitively on the initial 

state and the coupling to the environment. Entire decoherence- free subspaces (DFS) can exist if the coupling operators 

to the environment have degenerate eigenvalues [20l - |24| . 

We show below that a collective coupling that depends on a parameter x oi N quantum systems Si to a common 
"environment" TZ can be used to measure x with an uncertainty that scales as 1/N with an initial product state of all 
subsystems. The method works whether TZ is entirely under our control, or a reservoir with many degrees of freedom 
to which we have only partly access, i.e. a collective decoherence process of the Si, as long as we can measure an 
observable of the environment. 



Results 

Model 

Consider A'' quantum systems Si coupled to a common environment 7?.. The hamiltonian of the total system has 
the form 

JV 

H{x) - Yl ^' + Y. ^^-(^) ®Ru + hr, (1) 

■i— 1 i^u 

where Hi is the hamiltonian of system 5i, and for simplicity we take the Si as non-interacting. Hji denotes the 
hamiltonian of TZ, which may be itself a composite quantum system. Hamiltonian ([l} can be a model of decoherence 
(in which case TZ would be the ensemble of many degrees of freedom of a "reservoir" to which we have only partial 
access), or H(x) can generate a unitary evolution if TZ and S = {Si, . . . ,Sm} are completely under our control. The 
sum over v runs over an arbitrary number of operators for each subsystem Si and TZ, but R^, can also mean operators 
on different subsystems oiTZiiTZ is composite (e.g. positions of harmonic oscillators modelling a heat bath). In order 
to have a generic name for TZ that encompasses these different situations, we will refer to TZ as the "quantum bus" . 
The entire dependence on x is included in the coupling operators Si^i,{x). With "collective couplings" (and with 
"collective decoherence" if 7?, is a reservoir) we mean Si^i,{x) which do not depend on i. 

The smallest uncertainty 5x with which x can be measured is found from quantum parameter estimation theory 
[^. If the state of a system is given by a density matrix p{x), the smallest achievable 5x is given by 

6x > feinin = r-pr , (2) 

where we allow for M repetitions of the same measurement in identically prepared states p(x), and ds^ is a metric 
on the space of density operators. It is related to the Bures' metric dBmcsip, P + dp), with dp ~ p'[x)dx, by ds = 
2<iBurcs(/57 P + dp). For pure states, the Bures distance reduces essentially to their overlap, (iBmos(|'/')('0li I0)(</'I) = 
V2\/l — \{iIj\(J))\ [25|. If p{x) and p{x) + dp are related through a unitary transformation with generator h, p{x + dx) = 
exp{—ihdx) p{x) exp{ihdx) , then [2| 

dBnr:csipix),p{x)+dp) = {AhY^^dx. (3) 

An operational definition of 5x is given by 

(AA2)i/2 
6x = J,,,,/;,, , , (4) 



M\d{A).,/dx\ 

where we see that 6x corresponds to the quantum uncertainty of an observable A in state p{x), suitably translated by 
the slope of {A)^ into a fluctuation of x. As usual, for any observable A, (AA^) = {A"^) — (A)'^, and all expectation 
values are with respect to p{x). Inequality ^ holds for all possible measurements, and for A/ — > oo a measurement 
exists that saturates the bound Q. 

Model Q cannot be solved in all generality. However, if the interaction is sufficiently weak it can be treated in 
perturbation theory. Since we start in an initial product state at t = 0, we can then relate properties of the full model 
to the single particle dynamics of all Si and TZ. We first establish the fundamental lower bound on dx with the help 
of quantum parameter estimation theory, and then calculate dx for a given measurement on TZ. 

Quantum parameter estimation theory 

We decompose H{x) = Hq + Hi{x), Hj{x) = J^i v Si^^ix) (S) Ri^, and switch to the interaction picture with respect 
to Ho, with wave function \tpi{x,t)) = cxp{iHot)\'ip{x,t)), \'ip{x,t)) = exp(— ii7(x)i)|'0o)- In "Methods" we show that 
the Bures distance between the two states p{x) = \ipi{x,t)){'ipi{x,t)\ and p{x + dx) = \'ipi{x + dx,t)){ipi{x + dx,t)\ is 
given by 

dLres{pix).P{x + dx)) ^ dx" j [ dhdt2K\^,){H'j{x,h),H'j{x,t2)) (5) 

JO Jo 



where we have defined the correlation function for any two operators A, B in the state |V'), ^|i/j) {A, B) — {iIj\AB\'4)) — 
{il)\A\il)) {il)\B\^) , Hj(x, t) = exp{iHQt)Hj{x) exp(— iiJgi) is the interaction hamiltonian in the interaction picture, and 
Hj{x, t) = dHi{x, t)/dx. Equation ([5]) generahzes (|31), which is recovered if [H'{x), H{x)] — and Hj{x)t = xh. From 
([2]). wc have 
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(6) 



For identical and identically prepared systems Si, Si^y = S^, and |(p)i = |<^) for all i, and an initial product state 

l^o)-|v')^^®IO, (7) 



we find 



K\^,){H'j{xM).H'j{xM)) = E (NK\^^{Sl{xM).S'^[^M)){R.{tl)R^.{t2)) 



+N^Sl{xM)){S'{xM))K\^){Ru{ti),R^.{t2)) 



(8) 



where the expectation values for operators of Si (TV) are taken in states \ip) (|^)). Together with eq. ([6|) this proves 
the existence of a measurement on S and TZ that gives a 1/A^ scaling of 6x„^in ioi N ^ 1 and an initial product state, 
provided Eu. lo Io(SUX't^))(S'^(^'t^))K\oiRAti),R^{t2))dhdt2 ^ 0. 



Measuring the quantum bus 

We now use directly eq. ((4]) for showing that the 1/iV scaling can be achieved with the measurement of almost 
any observable A on TZ alone. The expectation values in eq. (j4]) are in general time-dependent. This implies a time- 
dependent minimal uncertainty as well which does, however, not affect the scaling with N. We evaluate {A{t)} and 
(Aj4^(t)) again by using second order perturbation theory in the interaction. The general results for these expressions 
are cumbersome, but simplify considerably if we make the following two assumptions: (1) The initial state of TZ is 
an eigenstatc of A, A\£^) ~ a^|^); and (2) A commutes with Hji. Both assumptions taken together imply that the 
quantum bus is prepared in a noiseless state at i = ((AA^(O)) = 0). Under the above two assumptions, we find 



(9) 

(10) 
(11) 



{A{t)) = a^+ dh dt2XS.^.{N,x,ti,t2)CRA.^.{ti,t2)+0{H'f), 
Jo Jo 

XSu^.{N, x,ti,t2) = NK\^) {S,{x, h), S^ix, ta)) + N^iS^^, h)) (5^(x, ^2)) 
CRAuf.{tut2) = {R,ih)AR^{t2)) - a^{RAti)Rt.{t2)) , 



which can also be used to obtain {A'^{t)) and (AA^(i)). Eq.(|3|) then leads to 
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IodtiJodt2j:.,^Xs.AN,x,h,t2){[R,{h),A][A,R^{t2)]) 



1/2 



JodhJ^dt2j:,^^£xs.^.iN,x,h,t2){R.ih)[A,R^it2)]) 



(12) 



In the limit of A^ 3> 1, the term quadratic in N in xsu^J.{N, x, ii, ^2) dominates, and we find a 1/A^ scaling of Sx, 

1/2 



5x 



Jo dh J^ dt2j:.JSAx,h)){S^ix,t2)){[RAti),A][A,R^it2)] 
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Jodhj;;dt2E.,^li{SA^,ti)){s,(.x,t2))){RAti)[A,R^it2)]) 
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(13) 



provided that the denominator does not vanish. It is enough to measure an observable of the quantum bus TZ alone, 
with all subsystems initially in a product state. 



Decoherence 

All derivations so far apply perfectly well if TZ is an environment with many degrees of freedom which wc cannot 
fully measure. Measuring an observable A on only a subset of these implies a non-unitary evolution of S. This 
establishes immediately that we can reach a 1/N scaling of 6x, if x parametrizes a collective decoherence process, and 
if we can measure at least some part of the environment. The example of superradiance that we will work out below 
is of this type. However, one might also be interested in how the unitary evolution generated by the hamiltonian ([T]) 
is affected by additional independent decoherence of the components Si and TZ. For Markovian decoherence, such a 
situation is described by a master equation for the density matrix W{t) of S and TZ of the form 

Wit) = -i (Lo + Li{x) + iKs + iAflO W{t), (14) 

where L^X ~ [Hq^X]^ Li{x)X = [Hi{x),X], and A^; (A5) are Liouvillians of the Lindblad-Kossakowski type (2a | 

for7^(5),withAs=EtlA^• 

The free evolution {Hj = 0) still factorizes, such that, essentially, all expectation values and correlation functions 
are replaced by expectation values with respect to the relevant mixed states (see eq. ([57)) in "Methods"). The 1/A^ 
scaling is therefore robust under individual decoherence of the components, an eventually increased prefactor not 
withstanding. This is corroborated by further exact results for a pure interaction with decoherence added to all St or 
to TZ (see "Supplementary Discussion"), and by the example of superradiance below. 

Measuring the length of a cavity 

As example of an application, we now show how to measure the relative change of length SL/L of a cavity with an 
uncertainty of order 1/N with an initial product state of N quantum resources. We first consider unitary evolution. 

Let N two-level atoms or ions (N even, ground and excited states |0)j, |l)j for atom i, i = 1, . . . , N) be localized 
in a cavity, and resonantly coupled with real coupling constants gi to a single e.m. mode of the cavity of frequency uj 

and annihilation operator a (see Fig.[T|), interaction hamiltonian Hj = 'J2i=i 9i (c- "^ + f^! a), where cr_ ~ |0)i(l|i, 

cr!|! ~ \l)i{Q\i. Due to the spatial dependence of the e.m. mode in resonance with the atoms, the gi depend on the 
position Zi of the atoms along the cavity axis and on the length L of the cavity (the waist of the mode is taken to be 
much larger than the size of the atomic ensemble) , 



5i = J— rFsin(/c^Zi)e-d, (15) 

where kz = 'kUz/ L^ eg denotes the dielectric constant of vacuum, V = LA the mode volume (with an effective 
cross-section A), e the polarization vector of the mode, and d the vector of electric dipolc transition matrix elements 
between the states |0)i and |l)i, taken identical for all atoms. 

If all gi are identical, we obtain the Tavis-Cummings model [27|. Here we consider the situation where the atoms 
can be grouped into two sets with N/2 atoms each and coupling constants Gi in the first set (i G {1 . . . , A^/2}), and 
G2 in the second set [i € {N /2 + 1, . . . , N}). One way of obtaining two coupling constants may be to trap the atoms 
in two two-dimensional lattices perpendicular to the cavity axis (see Fig. [1} . Note that it is not necessary to locate 
the atoms within a quarter wave-length of each other in order to obtain a DFS, as would be necessary without the 
cavity [2^. Distances which are integer multiples of the wave-length work just as well. In (jisp we have neglected 
the transversal dependence of the mode, assuming that the atoms are localized at a distance from the cavity axis 
much smaller than the waist of the mode. However, this is for a computational convenience only. The initial product 
DFS states also exist if there is a radial variation of the gi, but describing the dynamics would become much more 
complicated as it would depend on all the different gi. Assuming two different sets of coupling constants, the system 
is described by eq. ((!]), where we identify a pair of atoms (i, i -\- N/2) with subsystem 5^, i = 1, . . . , N/2 = Np, and 
the resonant cavity mode with the quantum bus TZ. The free hamiltonian Hq consists of the energy of all atoms. 
Hi ~ {uj/2){az + az ), and the energy of the cavity mode, Hn = ojo^a. 

An expansion of the gi about L for a small change 5L allows one to write the coupling in the form of ([1]) with 

Si,i{x) = g ({1 + x)a_ -I- (1 — x)a_ ) > '^'1,2 = S]-^, and Ri — a) , R2 ~ a, and x ex SL/L (see SI for the prefactor). 

For notational simplicity we restrict ourselves to the case where for x = the couplings are the same for the two sets, 
Gi = G2 = g, but this is by no means necessary for the method to work. 

A convenient basis for a pair of atoms is given by the "singlet" and "triplet" states {|s), |i_), |to), |i+)} with 
\s) = (|01) - |10))/y2, |i_) = |00), \to) = (|01) + |10))/V2, and \t+) = |11). As initial state of all 5, and TZ we take 




FIG. 1: Scheme for measuring the change of the length of an optical cavity. N atoms (or ions) are trapped at fixed positions 
in two 2D optical lattices perpendicular to the cavity axis. A dipole transition of the atoms is in resonance with a single, leaky 
cavity mode. The atoms are initially prepared in a dark state in which destructive interference prevents the photons from being 
transferred from the atoms to the cavity mode. When the cavity length L changes by a small amount 5L, the true dark states 
evolve, and the initial state is exposed to collective decoherence, detectable by photons leaking out through the semi-reflecting 
mirror at a rate proportional to A'^'^. This allows to measure 5L/L with a Heisenberg limited uncertainty of order 1/N, even if 
the initial dark state is a product state. 



the product state with N -^ N/2, and \(f) = (|i- 
We obtain a timc-indcpendcnt K{ti,t2) = g^ (Np 
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\s))/V2, and 10 = 
/2, and from © 



for a cavity mode in the vacuum state. 



6Xn 



MgW2N + N^ ' 



(16) 



vifhich clearly scales as 1/N ior N ^ 1. One might argue that a small amount of entanglement is present in \ip), but 
the size of the cluster of atoms all entangled vifith each other (i.e. a pair of atoms) is independent of N, such that it 
is legitimate to consider a pair of atoms as individual subsystem, and it is a product state of these subsystems that 
we consider. In SI we show that the product state can be prepared by letting the atoms interact pairwise. 

The initial state contains half a photon per atom. For a generic state the excitations stored in the atoms would 
start oscillating between the cavity mode and the atoms. However, for x = our initial state is a "dark state" , as 
destructive interference prevents the transfer of the photon from any pair of atoms to the cavity. When x deviates 
from zero, the perfect cancellation in the destructive interference is broken, and photons get transferred to the cavity. 

Measuring the number of photons constitutes an optimal measurement in the sense that the bound (|16|) is reached. 
To see this, we identify A = a' a in eq. ([T3| . This leads in a straightforward manner to Sx = ^^ , which agrees 
with ([TS]) for TV :^ 1, including the pref actor. After what was said in Sec. "Decoherence", it is clear that adding 
independent decoherence to all subsystems does not change the 1/N scaling of Jccmin. We now show this explicitly by 
considering the situation of very strong damping of the cavity mode, the superradiant regime. 

The framework of Sec. "Decoherence" is suited for this analysis, but we adopt the well-developed theory of super- 
radiance |29l - [3a | to give an independent demonstration that Sx^in scales as 1/A^. Decoherence arises because of two 
processes: Each atom can undergo spontaneous emission with rate F, due to its coupling to a continuum of additional 
e.m. modes. The damping of the cavity mode arises from the escape of photons with a rate 2k through one of the 
mirrors. In the notation of eg. lfHl). and identification of a pair of atoms {i,i + N/2) with Si, the generators A^ and 
Afl for these two processes read J29l - l32| 



ArX = K {[aX, a^] + h.c' 



(i+N/2) y ii+N/2) 
A,cr_,_ 



] + h.c) 



(17) 
(18) 



Superradiance occurs in the overdamped regime F ^ gvN <^ k, where a photon transferred to the cavity leaves the 
cavity before it can feed itself back to the atoms, but induces emission in other atoms while in the cavity mode. Cavity 
decay is then the by far dominant process. We will therefore start by neglecting F, but treat spontaneous emission in 
SI. The population of the cavity follows the occupation of the atoms adiabatically, and one can eliminate the cavity 
mode. This leads to the well-known and, for a; = 0, experimentally verified master equation of superradiance [23-[3J| 



for the reduced density matrix ps of the atoms in the interaction picture, 



-Ps{t) = Li{x)[ps{t)] = 7 {[J_{x)ps{t), J+(x)] + [J^{x),ps{t)J+{x)]) 



(19) 



The coUective generators J± arc J-{x) = Y.fli ^lAx) = X^^^/^^ ((1 + ^)'^-' + (1 - 2;)CTi'+^/^^V J+{x) = f_{x). 

The rate 7 = g^ / k, is independent of N . CoUective decoherence is a two-stage process here, as photons stored 
in the atoms first need to be transferred to the cavity mode before they can leave the system. The dark states 
of Sec. "Measuring the length of a cavity" are therefore decoherence-free states. There is a large decoherence-free 
subspace (DFS) containing (^2) ^ 2^ /^/N DF states, including a 2^/^ dimensional subspace ^iJi {\t-)i, \s)i} in 
which the pair formed by the atoms I and I + N/2 can be in a superposition of |t_); and \s)i [23, |35|. A DFS of the 
same dimension also exists for non-identical couplings, but the coefficients in the linear combination of the singlet 



state need to be adapted accordingly, \s) 



VigTF+Ig^ 



= (Gi|0);|l)i+jv/2 - G2|l)/|0);+7v/2)- If after preparing the 



atoms in a DFS state corresponding to the initial couplings G\ , G2 the length L of the cavity changes slightly, the 
coupling constants will evolve, G} — )■ Gj, I = 1,2, and so will the DFS. Photons will leak out of the cavity as the 
original state becomes exposed to decoherence. 

There is a well-known connection between the photon statistics in the cavity mode and the excitation of the atoms, 
derived in [301 fo^ all couplings identical. 



(a^™a"(t)) =2to(^) 



2m 



ds I 



'(e''''-l)2™-i(j:^J!"(t-s)) 



(20) 



One checks that this relation remains valid for small asymmetries a; 7^ 0. Thus, instead of the number of photons in the 
cavity for a given value x one can calculate the excitation of the atoms, where, however, the observable itself becomes 
a function of x, J™(a;) J™(a;). For the initial product state ([7|) considered above (with A^ -^ Np, \ip) = (|i-) + |s))/V2, 
and 1^) = |0)) we have from ^, 



(J+J_(i)) = 



'^ + ^ I - 7^ (3iV2 + A.3) (^2 ^ ^4) ^ 0(,2) ^ 



{jUUt)) 



2jt 




(21) 



-3Np + -Np 



0{t'). 



(22) 



Equation (|^n|) is in principle valid only in the Markovian regime t ^ 1/k, if {J^JZ''{t — s)) is obtained from 
the solution of the Markovian superradiance master equation ([T^. However, the initial behavior of (a^a(t)), 
{a^a{t)) ~ 5'^t^(J+ J_(0)), is entirely determined by the value of (J+ J_(t)) at i = 0, i.e. the question of the Markovian 
approximation of the dynamics of (J+J_(i)) does not arise, and eq. (|20p can therefore be used to calculate (nph(a;,i)) 
for short times up to order t^. From (J+J_(0)) one finds immediately 



(nph(.T,t)) = (a^aix^t)) = -gh^x^Np{Np + 1)) + Oit"") . 



(23) 



At this order k does not intervene yet, as initially the cavity mode is in the vacuum state. The quadratic initial 
increase of (^ph) reflects the beginning of a Rabi oscillation between the excited atoms and the cavity mode. We 
expect this result therefore to be valid as long as (nph(a;,t)) < 1. Equation (|23|) agrees identically with the result one 
finds from the approach in Sec. "Decoherence" (see eq. ([37| in "Methods"). The fluctuations of nph are obtained from 

(nph(a;,t)2) = {a'^^a'^{x,t)) + {nph{x,t)). Together with ((20]) one gets for {nph{x,t)) < 1, (An^j^(x,i)) ~ {nph{x,t)) 

with corrections of order (q/n)^. From eqs.dH) and (1^51) we flnd Sx ~ ^= — ^^ ~ jC^ — , which is identical to 

the minimal possible uncertainty, ea. p6)) for A^ ^ 1. The validity of the short time expansions pil22p is limited to 
A'^7f <C 1, as can be seen from comparing the first order term with the zeroth order term. Inserting (|2ip in (j20p gives 
therefore an analytical prediction of (nph(O) valid for g/n <^ gt <^ K/{Ng), in addition to the small time result (|25)) 
for gt <C g/n- The agreement of (nph(i)) based on (|21l) with the result from simulating (jl9|) can be further improved 



by rc-cxponcntiating {J+J-{t)) according to a + ht ~ acxp(6/at), before inserting it in (jSO]) . The limitation of validity 
of the small-time expansion does not pose a serious restriction in the bad cavity limit k ^ g, nor does it imply that 
the 1/A^ scaling of the sensitivity breaks down beyond that regime. A full theoretical analysis for longer times will 
have to include the calculation of the superradiant propagator with broken SU(2) symmetry, however. For t ~ l/n a 
non-Markovian description of superradiance is called for, which is beyond the scope of the present investigation. 
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FIG. 2: Mean photon number and uncertainty Sx of the change of length of the cavity (a): Mean photon number (riph) as a 
function of dimensionless time gt (where g is the couphng constant of the atoms to the cavity mode) for A^ = 2, 4, 6, 8, 10, 12 
(black, red, green, blue, brown, violet), in units of x^ for x — 0.1, and with photon escape rate from cavity k = 5g, obtained 
through numerical simulation of superradiance using a stochastic Schrodinger equation. The dashed lines with corresponding 
colors are analytical results valid up to Ng^t/K ~ 1 (see eqs. (|20p . (|2ip V (b): Uncertainty 5x (see eq. Q) based on A — Tiph 
as function of A^ for gt — 0.0485 and x — 0.01. Numerical results (circles) show the same 1/A*' scaling as the ideal lower bound 
(red dashed line), eq. p6[) . with slightly increased prefactor. Green continuous line is an analytical prediction based on an 
expansion of {nphit)} for small gt. 



Figure [H shows that (nph(t)) obtained numerically by simulating (|19l) through an equivalent stochastic Schrodinger 
equation (SSE), and integration of (J+J_(f)) according to eq. (l20l) . agrees well with the result based on ((2T|) for 
gt < g/K and g/n < gt < K/{Ng). The SSE for real ip(t) reads 



diit) = Di{iit)) dt + D2{iit)) dW{t) , 



(24) 



with Di(V^) = 7 ('2(J_)^ J_ - J+J^ - {J-)'^] 'ip, and D2(V') = V^ {J- - {J-)i^) ip, where dW{t) is a Wiener process 

with average zero and variance dt, and {J-),p = {ip\J-\ip) [2g. We used 2000 equidistant time steps in the time 
interval i = 0, . . . , 20/ g, 20 random realizations of the process for the simulation of (nph(i)), and 400 realizations for 
the calculation of Sx. Figure [2] also shows 5x calculated from the numerical data for (Anpj^(i))^/^ and {npi^{x,t)) 
through eq.(|l]), together with the fundamental lower bound Sx^in, ea. ([T5|) . We sec that at gt = 0.0485, 6x follows the 
optimal 1/A^ scaling with only slightly increased prefactor. 

We emphasize that Uph allows to measure SL/L, not just to detect a change of L. Egs. ipO)) and (PT|) relate (nph) to 
X, and, unless the two lattices are situated at anti- nodes of the mode, the relation between SG = (Gi — G'2)/2 = gx and 
6L/L is linear to lowest order and independent of N: If we choose the position of the atoms such that Z2 — zi = m\ 
with riz — 1 > m G N we have 



SG /rirTTZi 

— = mn cot -^ 

g \ L 



SL 



Therefore, the measurement of (riph) allows the measurement of SL/L. Several other practical questions, e.g. the 
preparation of the initial state, and the robustness of the method with respect to fluctuations of the coupling constants, 
spontaneous emission, and errors in the preparation of the initial state, are addressed in SI. The superradiant regime 
has the advantage of providing direct access to the number of photons in the cavity. The average number of photons 
outside is simply obtained by integratin g 2K {a^ a{t)) up to time t, as the photon escape rate is proportional to the 
average photon number inside the cavity |30|. The results for the scaling of Sx with N based on ^4 = a^a are therefore 



unaffected by detecting the photons that leave the cavity. Measuring the number of photons amounts to monitoring 
the decoherence dynamics, and we have thus an example where the parametric dependence of a collective decoherence 
process allows to achieve the Heisenberg limit with an initial product state. 

Discussion 

Our results may seem to conflict with the well-known theorem [J] that for unitary evolution of N independent 
quantum systems in an initial product state at best a scaling (5xniin ^ ^/vN is possible. To see that there is no 
contradiction, it is helpful to consider the simple case where H'{x) = dH{x)/dx and H{x) commute, [H [x) , H' [x)] = 0. 
One then easily shows that to lowest order in dx, p{x, t) = exp(— iiJ(x)t)p(0) exp(iiJ(a;)t) (with fi, = 1) and p(a;, t)+dp 
are related by a unitary transformation with generator h = H' {x) t. Let us furthermore restrict ourselves to a single 
operator per subsystem, i.e. v = 1 only, and to the linear x dependence 5*^^1(2;) = xS for all i, and Ri ~ R. A few 
lines of calculation lead to 

(Ah^) = {N{AS^){R^) + N^{S)^AR^)) f , (25) 

for an initial product state, p(0) = |V'o)('0o| with jV-'o) from eq.©. All expectation values of S in (|25|) are in state |(p), 
those of R in the state |^). Inserting ^ into © and ©, we find that for iV > 1 and {S)'^{AR^) ^ 0, 



2VMi|(S')|(A,i?2)i/2iV 

i.e. the Heisenberg limit 5x^i^ ^ 1/N can be achieved with an initial product state. Clearly, for the case considered 
above the unitary transformation generated by Hq = ^,- Hi + Hr is not necessary in order to achieve the 1/N scaling. 
We therefore simplify the reasoning further by considering the case Hq ~ 0. We are then left with a pure interaction, 

H{x)^Hi{x)^xJ2Si®R- (26) 

i 

But this is not a hamiltonian of the form H(x) = x ^^ hi required by the theorem in (j]. In our case all subsystems 
couple in a non-trivial fashion to the common quantum bus TZ and are therefore not independent. This turns out to be 
the decisive difference. The SQL can be recovered for the standard situation of N independent subsystems through a 
R that acts only trivially on 7^, i.e. i? = 1, such that {AR^) = 0, and thus fe^in = 1/ {2y/ M Nt{AS'^)^/'^). This makes 
obvious the rather ironic fact that quantum fluctuations in TZ help and are necessary to achieve the 1/A^ scaling. The 
prefactor of the 1/A^ behavior is smallest for an initial state with an equal weight superposition of the eigenstates of 
R pertaining to its largest and lowest eigenvalues rmin and rmax, in which case (Ai?^)^'^ = l^^max — ''min|/2. 

This simple example also allows to corroborate that in order to achieve the 1/A^ scaling one need not measure the 
Si at all, and almost any measurement on TZ suffices. Consider an initial product state l-^o) = '^iLi\si)i\Oj with 
Si\si) = Si\si), and \£,) = J^md-mlrm) for i?|r™) == r™|r„i). We then have 

N 
m i—1 

Let A be an observable on TZ which does not commute with R, i.e. there are at least two eigenstates |ro) and \ri) 
such that (ro|A|ri) 7^ 0. It is sufficient to consider an initial state of TZ which is a superposition of these two states, 
e.g. we may take do — di = l/"\/2, and an observable A ~ |7'o)(ri| + |ri)(ro|. If all subsystems Si are prepared in the 
same state with Si = s, one finds {A{t)) = cos(a;iVs(ro — ri)t) and {AA^{t)) = sn:?{xN sir^ — ri)t). Inserted in eq. ^ 
this leads to the exact result 

A^|s||ro - ri\t 

valid for all x. This shows that the 1/A^ scaling can be reached by measuring almost any observable of 7?., as long as 
it does not commute with R. Furthermore, eq. (|27p allows a simple quantum information theoretical explanation of 
the effect: The final state reflects the accumulated phase from the interaction of all the systems Si with the common 
quantum bus TZ. Figure [3] shows an equivalent quantum circuit that reproduces state (|27p . One subsystem Si after 
another imprints the same phase on the components of the state of TZ. Equation (j27p also makes obvious that a 
measurement of the Si alone does not allow to achieve the 1/N scaling, as the state of TZ will collapse on a single 






ro) + \ri) - H,,i(x) - H,,2(x) - H,,3(x) H,,a/(x) - |ro) + e'-" |ri) 



FIG. 3: Quantum circuit that reproduces the interaction hamiltonian Hi{x) — x '^^^^ Si®R. N quantum systems Si prepared 



all in the eigenstate \s) of Si, Si\s) = s\s) in eq. (|26p . lead to a total accumulated relative phase between states of the quantum 
bus 7?. that is proportional to xN. This allows a measurement of x with a precision that scales as 1/A*', even though the initial 
state is a product state, by measuring any observable A on 7?. alone that does not commute with R. 



state \rm), and one only gets an irrelevant global phase. Thus, measuring the quantum bus is not only sufficient, but 
also necessary for the 1/A^ scaling of 6x. We also see that the measurement oi A = \ro){ri\ + |ri)(ro| is optimal if 
ro,ri correspond to the smallest and largest eigenvalues of R, respectively. 

In |13j an adaptive measurement technique was demonstrated that allows one to achieve Heisenberg-limited uncer- 
tainty by using only an initial product state. The method is based on phase estimation |36| . but instead of using a 
NOON state of N photons, independent photons were passed N times through the same phase shifter. This amplifies 
the phase by a factor TV, but it was shown that in the presence of losses the scaling of the sensitivity with N is at 
most improved by a constant factor 37| compared to the classical case for TV — > oo. A common feature of both phase 
estimation and our method is that a measurement is performed on a common quantum system that interacts with 
all other quantum systems. However, our method is more general. It incorporates decoherent and unitary evolutions 
in the same framework, and allows one to use collective decoherence as a signal. Secondly, phase estimation was 
developed for a multi-qubit system with controlled, sequentially turned on interactions, and an a;-dependence in the 
free evolution. Hamiltonian ([T]) on the other hand can be used to describe substantially more complex systems, with 
possibly non-trivial dynamics in the absence of the collective interaction, and with interactions that do not com- 
mute with the Hamiltonians of the free constituents. Furthermore, the interaction is simultaneous such there is no 
bandwidth penalty in the accumulation of the phase, nor is there a need to re-sample a phase shift many times. 

Our method requires a collective interaction between A^ separable quantum systems and a common quantum bus, 
an initial noise-less state in the sense discussed above, and the possibility to measure at least part of the quantum bus. 
In the case of incomplete measurement of the quantum bus this implies the need of a collective decoherence process 
with a decoherence-free initial state. Besides atoms in a cavity one might consider circuit-QED systems [Sg , trapped 
ions coupled to a common phonon mode [3^ , or quantum dots coupled to micro- resonators [40| or to photonic crystals 
[41| . Both unitary evolution or a decoherence process can be useful, as long as the collective interaction between the 
N quantum resources and a common quantum bus depends on the parameter x to be measured. 

To summarise, we have developed a general theory of collectively enhanced quantum measurements based on the 
interaction of N quantum systems with a common "quantum bus" . The latter can be a simple quantum system, or an 
environment with many degrees of freedom to which we have only partial access. We have shown that if the collective 
interactions depend on a parameter x, the Heisenberg limit (i.e. a 1/iV scaling of the uncertainty of x) can be reached 
with an initial product state, and by measuring almost any observable of the quantum bus. We have used quantum 
parameter estimation theory to establish that a 1/7V scaling of the uncertainty is indeed optimal in this setup. We 
have given a simple quantum- information theoretical interpretation of the effect, and we have analysed in detail a 
possible experimental implementation of the measurement of the change of the length of a cavity with an uncertainty 
that scales as 1/A^. 

The proposed measurement principle offers an attractive way out of the dilemma of ubiquitous decoherence that has 
so far plagued quantum enhanced measurements: First of all, there is no need to build highly entangled states which 
are extremely fragile under decoherence for large N. Simple product states will do, and decoherence of some parts 
of the system does not affect the 1/iV scaling of the minimal uncertainty. Secondly, parameter-dependent collective 
decoherence is covered itself by our new measurement principle. Indeed, decoherence is a process in which quantum 
interference effects can play an important role. This is exemplified by the very existence of DFS, and can lead to 
exquisite sensitivity when a DFS is disturbed. Instead of trying to suppress decoherence at all costs, one might 
therefore be better off exploiting its parametric dependence. 
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Methods 

Bures distance for unitary evolution 

The state vector |V'/(a;,t)) in the interaction picture obeys the time-dependent Schrodinger equation 

i— |V'7(x,t)) = Hj{x,t)\yjj{x,t)) , (29) 

with the interaction hamiltonian 

Hi{x,t) = J2^iAx,t)®R^it), (30) 

S,,,{x,t) = e'^-*5,,,(a;)e-'^-*, i?,(t) = e'^«*i?,e-'^^* . (31) 

The general solution of ([29]) is given by \ipj{x,t)) = Tcxp — i/p Hj{x,t')dt' |'!/'o)j where T denotes the time-ordering 

operator. To second order in the perturbation Hj{x,t), the overlap between \ipj(x,t)) and \ipi{x + dx,t)) reads 

pt pt pt\ 

{iji{x + dx,t)\iJi{x,t)) = 1 + dxi {H'j{x,ti))dh+ / {[H'j{x,ti),Hi{x,t2)])dtidt2 

L Jo JO Jo 

rt pt pti 



ldx^(i I {Hj{x,ti))dh+ I I {[H'j{x,ti),Hi{x,t2)])dhdt2 
^ ^ Jo Jo Jo 

2/ r {H'j{x,h)H'j{x,t2))dhdt2) , (32) 



Jo 



with all expectation values with respect to \ipo). We assume that the derivatives of Hj{x) with respect to x are 
hermitian operators, in which case the term linear in dx is purely imaginary. The lowest order term in the squared 
overlap is then of order dx^ . One finds in a straightforward manner the squared Bures distance ([S]). 

Decoherence of subsystems 

Markovian decoherence of the Si and TZ on top of the unitary evolution generated by H{x) can be described 
by eq. (HH). The single system dynamics (i.e. all Si and TZ taken separately, Lj = 0), can be solved formally by 
exponentiating the Liouvillians. Wc will again treat Lj{x) in perturbation theory. The density matrix Wi{t) of S 
and TZ in the interaction picture is related to the one in the Schrodinger picture, W(t), by 

Wit) = e-'(^o+'A=+'A'-)*VK/(f) = e-'^^^'Wiit) , (33) 

and obeys the master equation Wi{t) — —iLi{x,t)Wi{t). With eq. ([55]) we have defined the free propagator Ppit) = 
e"'^* = (8)ie~'^'* (g) e~'^«* = ®iPi{t) ® Pnit), and Lj{x,t) = PF{-t)Li{x)PF{t) is the interaction Liouvillian in the 
interaction picture. We decompose furthermore Li ~ X]fc=i ^i.ki with Lj^X = [Hi^k,X]. To second order in i/ we 
have 

{A{t)) = tr{A(Ff(i)-i^ / Pj.(i-ti)L/,fc(x)Pf(ii)dii (34) 

- y, f dh r dt2PF{t-h)Lj,k{x)PF{tl-t2)Lj,k'{x)PF{t2)+0{L]t,))W{Q)]. 

k,k'-'o J^ I i 

With an initial product state, W^(0) = (8'^i/5i(0)®/9r(0), we obtain the zeroth order term (A(t))o = \.T{APF(t)W{^)) = 
tr((8)iPi(i)[jOj(0)] (8) APR(t)[pfl'(0)]) = \.Yu{Apii(t)), as all propagators are trace-preserving, and puif) = Pr (i) [pflXO)] . 
Similarly, by explicitly writing LjkX = J2i,[SkMRi^,X], we obtain the first order term 

(A(t))i = -i^Ztr^A / PF{t~h)Lj^kPF{ti)dhW{Q)] (35) 

= -iV [ dti{Sk.Ati))trR{APRit - h)[R,, pnih)]) , (36) 

k,u-'° 
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where {Sk,,^{t)) = tikSk,vPk{t)[pk{0)]- We generalise the second simphfying assumption in Sec. "Measuring the 
quantum bus" to PR{t)[pR{0)] = pr{0), and tr R{APR{t)[X]) = fA{t)tTR{AX) with some function /^(i) H^- This 
imphes that the initial state is decoherence-free concerning the decoherence of Ti. alone. This is a natural assumption 
for the state of an environment initially in thermal equilibrium, or for a quantum bus in its ground state, such as an 
initially empty cavity mode (see the example of superradiance). We then have again {A(t))i = 0. To second order in 
the interaction we find 

(Ait)) = (A(0))o 

rt rt 

E 



u,t,--^ -^0 



+ N^{S^{t,)){S,{h))cf^^,{tMM)}, (37) 

CAl^^uit^M) = tTR{APRit^h)[R^,PRih-t2)[R,PRit2m (38) 

'^iflM-(*'^l'^2) = tTRiAPRit-h)[R^,PRih-t2)[pRit2)R.]]) (39) 

c'^l.^.At,ti,t2) = c^X^^itMM-cfR^StMM) (40) 

C^sl^tiM) = tr,. {S^Pk{t, - t2)[S,Pk{t2)]) - {S,ih)){SAt2)) (41) 

4'i(^l'^2) = tl-kiS^Pk{h-t2)[pkit2)S.])~{S,ih)){S,it2)), (42) 

where pR(t) = PR{t)[pR(0)]. The index k is arbitrary, fc = 1, . . . , A^, as we have assumed all systems Sk identical and 
identically prepared. All x dependence is in the operators S^. Equation ([57)) also gives (A^(i)) by replacing A ^ A^, 
and {A{t))^ to order 0{Hj). The equation obtained by inserting these expressions into eq. Qi generalizes the result 
([T2|) to decoherence on top of the unitary evolution considered in Sec. "Measuring the quantum bus" . We see that the 
basic structure of the result for Sx^ and in particular its scaling with N is unchanged, but the expectation values and 
correlation functions arc replaced by more complicated expressions involving in general mixed states and non-unitary 
evolution of individual subsystems. 

Quantum parameter estimation for a Markovian master equation 

In standard descriptions of decoherence one traces out the heat bath and gets a master equation for the reduced 
density matrix ps of 5 alone. Using quantum parameter estimation theory generalized to non- unitary evolution we now 
show that for Markovian decoherence with an initially decoherence-free state, measuring an arbitrary x-independent 
observable on S alone gives at best a Sxmin "^ ^/VN. This corroborates the result found for unitary evolution that 
the important quantum system to measure is the common quantum bus TZ, rather than S. 

The Markovian master equation for ps (t) obtained by tracing out TZ has the Lindblad-Kossakowski form 

d 

Ps{t) = A(a;)[p,(t)] ^7^ ([i^„(x),p,(t)F2(a;)] +/I.C.) , (43) 

Q = l 

where we work in the interaction picture and assume that there is no additional unitary evolution. The Fa{x) are 
arbitrary linear (not necessarily hermitian) operators which have inherited the x-dependence from the interaction 
hamiltonian Hj{x), and d is the total number of generators. Note that we can restrict ourselves to an initially pure 
state, as for any linear propagation one cannot do better with a mixed state than with the pure states from which 
it is mixed ^^. We expand the Markovian time evolution to first order in t, p = |'i/')('0| + t ■^{x)\i') {il'l + 0{t^), 
and linearise Fa{x) about the value of x where we want to measure. We set that value, without restriction of 
generality, to zero, i.e. Fa{x) = xF^ + Fq,(0), and assume that the initial state is decoherence-free at a; = 0. The 
Bures distance can still be evaluated in a straight-forward fashion as the state at a; = remains pure. One finds 
d'^^ \x=o ~ ^T* Sa ^W i-^ai Pa) dx^ . As a conscqucncc, the ultimate quantum limit of the sensitivity with which the 
parameter x can be estimated from the parametric dependence of the master equation, starting from a pure state |-0), 
reads 

SXn,in = : ■ ^72 ■ (44) 



2V2M^ ELi^lv>>(^<i:^") 
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With Fa = Y.r=i^a.,r we obtain K\^){Fl,Fa) = Z]r,s=i -^l»A>(^lr'^a,s)- For an initially entangled state, 
^\ip) (-^ai ^a) can be of order A^^. This can be seen from the example of the GHZ state |V') = (|0 . . . 0) + |1 . . . l))/-\/2, 
and a single generator Fi^i ~ CTz with cri the Pauli z matrix for subsystem i. Then K\^^{Fl^Fi) — N^, and 
one obtains a 1/N scaling of (Jxmin, just as in the case of unitary evolution. However, if the initial state factorizes, 
\ip) = ^iLiI'pOj there are no correlations between different subsystems r and s, and we thus have only the sum of 
correlations in all subsystems, i^|^) {Fl,Fa) = X]s=i ^lv»> (^a si -fa.s)) which is at most of order N, and Sx^in scales as 
1/VN, again just as in the case of unitary evolution. This shows once more that a measurement of an x-independent 
observable on S does not allow to do better than in the standard situation of unitary evolution of S without coupling 
to a common quantum bus. 

Interestingly, superradiance is described by a master equation of S alone after tracing out the cavity mode. But a 
measurement on TZ (the number of photons in the cavity) translates in that case to a measurement on the Si that 
depends itself on x. In this way it is still possible to achieve a 1/A^ scaling of Sx. 

The authors declare to have no competing financial interests. 
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Supplementary Discussion for the manuscript "Heisenberg-limited sensitivity with 

decoherence-enhanced measurements" 



Pure interaction and decoherence 

We reconsider here the simple interaction-dominated model (26), but with additional decoherence of the subsystems. 
To simplify things and make connection to the language of quantum information theory, we consider directly the 
case where both the Si and TZ arc qubits. In the eigenbasis of Si and i?, the interaction Hj(x) can be written 
as Hj{x) = xJ2i=i ^z <8i fi , where we use the label for TZ. As observable on TZ we chose A = a^- Without 
decoherence one obtains the uncertainty Sx = l/{2\/MNt) in agreement with (28). We will restrict ourselves to 
phase- flip channels acting independently on all qubits with a rate F, 

A^„X ^ ^ ([a«A,aW] + [a«, AaW]) . (SI) 

As before we decompose the Liouvillian Lj{x)X = [Hj{x),X] ~ J2i=i^i,i'^- C^ne easily verifies [A (i),i/,fc] = for 
all i, k. 

Consider first phase flips of the qubits Si. The equation of motion for the full density matrix W{t), W{t) = 
'^i^iiLi,i + A (i))W, is easily solved in the computational basis, (a = a^aN^i . . . ao, ai = ±1, i = 0, ...,A^), 

Wah (t) = Wab (0) exp ( — i X] j=i (^ (q^jOo — bibo) + ir{aibi — 1)) i ) ■ The reduced density matrix for TZ alone is therefore 



5aofco = X! W^a„...aiaoa„...ai6o(0)exp I -i^a;aj(ao -5o)t j 

ai,...,aiv— ±1 \ i—1 / 



(S2) 



We see that T drops out completely, and dephasing of the Si has therefore no impact at all on the total accumulated 
phase of TZ, and therefore on Sx. 

If phase flips occur in TZ, W = [{J2i=i ^i,i) '^ ^cr^°n^^ ^^ have Wah{t) = 

I^ab(O) exp (-i (j2Zi ^ (a^ao - b,ba) + ir(aofoo - 1)) i) • For the initial state IV-o) = I + • • • +>(!+) + |-))/\/2, 

we find the reduced density matrix element p-| it) = exp(— 2rt) cos(2iVa:;t). Thus, the decay of {ax{t)) is independent 

of TV, and only the prefactor in the scaling of 6x with N changes. 

Superradiance: practical issues 

We now address several additional technical questions that are important for a possible experimental implementation 
of our method of measuring the length of a cavity using decoherence-enhanced measurements. 

Preparation of initial state 

In order to prepare the product state |^o) = {'^iLi-7^i\t~)i + \s)i)) |0) it is helpful to use three-level atoms with 

a lambda structure. Let |0) and the additional state |2) be hyperfine (HF) states, and assume that their energies are 
sufficiently split such that only the transition |0) <-^ |1) resonates with the cavity mode. We assume further that the 
second optical lattice can be moved along the cavity axis, such that controlled pairwise collisions of corresponding 
atoms in the two lattices can be induced. Entangled pairs of atoms in their HF split ground states can thus be created 
(for atoms in the same lattice this has been demonstrated experimentally, see Supplementary Reference [43] for a 
review). After the creation of an entangled HF state \iP'q), that differs from |-0o) by the replacement of states |1) by 
states 1 2), the second lattice is moved back to its original position. Now one can selectively excite the |2) states by 
a laser pulse in resonance with the |2) o |1) transition, that replaces the singlets in the (very long lived) HF states 
by the desired singlets of the |0) and |1) states and thus produce |i/)o). However, as such, the method is not of much 
practical use yet, as it will be virtually impossible to park the second lattice at the exact position corresponding to 
coupling constants which render |'0o) decoherence free. The extreme sensitivity of the collective decoherence with 
respect to changes of the coupling constants plays against us here, and will lead to leaking of light from the cavity 
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after the excitation jf/'o) -^ \ipo), if the exact position corresponding to j-iAo) € DFS is not achieved. But it is possible 
to position the second lattice at the required position with a precision of 0{1/N) using a feed-back mechanism and 
a part of the quantum resources. With the atoms in the state l^g), do the following repeatedly in order to find 
the optimal position: Excite a part of the entangled HF pairs containing 0{N/hiN) atoms with the laser, measure 
(riph), and use the measurement results to bracket the minimum of (nph) as function of the lattice position. The 
minimum of (riph) indicates that the position corresponding to the DFS is achieved. Using golden section search, the 
minimum can be bracketed to precision 1/N in C(ln JV) moves, as at each step the sensitivity of the measurement of 
the position of the lattice is of order 0{lnN/N) ~ 0(1 /TV). Once the minimum is found, excite the remaining unused 
pairs (there should be still a number of pairs of 0{N)) to the desired state \tpo)- That state is now decoherence-free, 
and the system ready to detect small changes of the position of one of the mirrors. Note that for this method it is 
not necessary to know which exact state is produced in the controlled collisions and subsequent laser excitation. 

Imperfections 

Here we discuss how the most important additional noise sources affect the photon statistics. In view of eqs. (20,21), 

we restrict ourselves to analyzing a ~ —{Jz{x,0)) = tT{JzLi{x)[ps{0)]) ~ 27(J+ J_(a;, 0)), where Jz = ^J2i=i'^z 
is the total population inversion of the atoms. Energy conservation shows that a can be interpreted as the initial 
photon escape rate from the cavity. 

Spontaneous emission 

It is easily verified that the spontaneous emission term in eq. (17) leads to a contribution to (7^(0)) that scales 
as 0{N), to be compared to the term of O(iV^) from collective emissions, cq. (21). The modification of the initial 
state due to spontaneous emission is not an issue, as we consider an initial product state. Atoms which decay become 
simply unavailable for collective decoherencc, but since their number is proportional to N this does not change the 
scaling of 6x with N. Also, note that spontaneous emission sends photons into the entire open space but not into the 
cavity, whereas the collective emission escapes exclusively through the leaky cavity mirror. Therefore, in addition, 
the two contributions can be well separated experimentally by observing only the photons which escape through the 
cavity mirror. 

Fluctuating coupling constants 

Another obvious concern are fluctuations of the coupling constants. In order to reduce the noise of the measurement 
of (wph), one may want to repeat the experiment M times with M ^ 1. The exact coupling constants might fluctuate 
during the averaging, e.g. due to fluctuating traps caused by vibrations in the set up. But even for perfectly stable 
traps, thermal motion, or even quantum fluctuations in the traps will lead to fluctuating g^. We now show that the 
cost in sensitivity of these fluctuations depends on their correlations. Fluctuations correlated between pairs of atoms 
/ and I + N/2 come at no cost, completely uncorrclated fluctuations lead back to a 5a; ex l/\/iV, and fluctuations 
perfectly correlated within the same lattice but not between the two lattices lead to noise indistinguishable from the 
signal. 

To see this, consider fluctuations 6gi of the gi about their mean values G/, gt = Gi + Sgt for i = l,...,A^/2, 

gi ^ G2+ Sgi for i = N/2 + 1, . . . ,N. The generator J_ reads then J^ = - X]i=i 9i'^- ■ We introduce the correlation 
matrix Cij = SgiSgj/g^, where the over-line denotes an average over the ensemble describing the fluctuations, and 
assume % = for simplicity. Then J-|V^o) = -^ Z^ili (5j - 9i+f)\t~)i (8>i^i \fi)i ^nd 



N/2 N/2 

2^2 



( ^(di - 9t+N/2)'^ ^ ^ 2 ^■^'' ^ •^'+^/2) (^J " 9J+N/2) J ■ (S3) 






Equation (jS3p predicts a background abg in the photon escape rate on top of the value for Sgi ~ 0, 

ao = ^x\2N + N^) , (S4) 
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with X ~ (Gi — G2)/(2g), sec eq. (21), and fluctuations Saf, a ^ ao + a^g + ^ctf, where the average background is 
given by 

JV/2 N/2 

1=1 i^j 



where Cij = ICij + C^^n_j+k — C^j^n. — Cj_|_iv j. The average background abg can be determined independently 

at X = 0, and subtracted from the signal; it does therefore not influence the sensitivity of the measurement. The 
remaining noise Saf fluctuates about zero, 

^«/ = x,(^+i)x:^, (S5) 

1 /9 

where Agi = Sgt — Sgi^i^/2- Its standard deviation {Sa'iy^'^ = 7|a;|I?; where D = (^- + 1) ( ^^ ^^-^ Cij j translates 
into additional noise in the number of detected photons. Several interesting cases can be considered: 

1. Fully uncorrelated fluctuations, Cij = CiSij, where Sij stands for the Kronecker-delta: Here we get D = 

1 /9 

(-J + 1) {J2i=i i^i + C'i+jv/2) ) , which is in general of order iV"^/^, and leads back to the standard quantum 



limit, Sx ^ 1/viV for large N. 

2. Pairwise identical fluctuations between the two sets: Cij = Ci^N_ j = C^jj^n_ = C^_^Nj_^n_ for i, j = 1, . . . , N/2. 
This can be the consequence of fully correlated fluctuations, Cij — C V i, j . Alternatively, such a situation arises 

for example for atoms initially arranged symmetrically with respect to an anti-node such that G\ — G2 , ^i 
the two atoms (or ions) in each pair I [l — 1, . . . , N/2) arc locked into a common oscillation. This should be 
the case for two trapped ions repelling each other through a strong Coulomb interaction, and cooled below 
the temperature corresponding to the frequency of the breathing mode. Equation (|S5[) then gives 5af = 0, 

i.e. no additional noise from the fluctuations of the couplings. Note, however, that for initial G]^ 7^ G2 the 
more general DFS leads to a more complicated condition for the correlations, GijjGj P + Gj^n aj^n_\G{ ^ — 
Cijj^N_G2 Gi — Cj^jv ^G-^ G2 — 0, which might be harder to achieve. 

3. Correlated fluctuations within a set, but uncorrelated between the two sets, Cij ~ C for i, 7 e {1, . . . ,N/2} or 
i,j e {A^/2 + l,...,A^},but C^j = Ofori e {l,...,N/2) and j S {N/2 + 1, ... ,N) or vice versa. This case leads 
to a noise of order 0{N'^), the worst case scenario. However, this comes as no surprise, as such correlations are 
indistinguishable from the signal: all the atoms in a given set move in a correlated fashion, but independently 
from the atoms of the other set. This modifies the couplings in the same fashion as if the length of the cavity 
was changed. 

Case (2) above is clearly the most favorable situation. If there are no other background signals depending on N, we 
keep the 1/N scaling of Sx. In order to favor case (2) over cases (1),(3), it appears to be advantageous to work with 
ions and to try to bring the ions in a pair as closely together as possible, thus strongly correlating their fiuctuations, 
while separating the ions in the same set as far as possible. 

Imperfections in initial state preparation 
Suppose that instead of the state \ijja), the state of the atoms 

N/2 
1=1 

was prepared (we consider the same state for all pairs for simplicity, but this is not essential, and assume the state 
normalized). Then, 

2 
-y (y - 1) ((G2 - Gi)&(a - d) + (G2 + Gi)c(a + d))'] /g^ . (S7) 



a = 7 



(((c - b)Gi + {c + 6)G2)' + 2d^{Gl + G^)) 
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The derivative of a with respect to Gi is of order 0{N'^), and thus still allows to find the minimum of a as function 
of the position of the second lattice with a precision of order 0{1/N). At the minimum a component outside the 
DFS persists, such that photons will leak out of the cavity, but the average rate is only of order 0{N), and can be 
measured separately and subtracted from the signal. Changes 5G2 of G2 away from the position of the minimum still 
lead to a signal that scales, for large N as N'^, a = j {a{c + b) + d{c — h)) N'^5G\/ g^, and the analysis leading to the 
1/A^ scaling of 5x still applies. 

Alternatively, one can get rid of the additional background by letting the system relax before measuring changes 
of L. Indeed, any state with a component outside the DFS will relax to a DFS state or mixtures of DFS states 
within a time of order I/7 or less. Components with large total pseudo-angular momentum J relax in fact in much 
shorter time of order 1/{J^). The DFS states reached through relaxation starting from ji/'o) still allow a scaling of 
a close to N"^ . We have shown this by simulating the relaxation process with the help of the stochastic Schrodinger 
equation (24). Using an Euler scheme with a time step of O.OI/7, we followed the convergence of V'(i) to DFS states 
for states with (a, 6, c, d) = (cos 5, 1, 0, sin(5)/-\/2, until the norm of the difference \^{t + dt)) — |V'(0) dropped below 
10""'^^. In these final dark states, randomly distributed over the DFS, we calculated a, and averaged over a large 
number n^ of realizations of the stochastic process (n^ = 10^, 10^, 1 0^, 2.5 • 10^, 10^, 1.25 • 10^, 250, 200, and 
250 for iV = 2, 4, 6, 8, 10, 12, 14, 16, and 18). [Supplementary Figure SI shows the scaling of a as function of TV for 
different values of 5 for < 5 < 'n 12 up to A^ = 18. Within this numerically accessible range of A, a follows a 
power law a ex A^^ with an exponent p that decays only gradually with b for b < 7r/4. Moreover, that decay might 
be a finite size effect: Note that, surprisingly, a{S) appears to be close to symmetric with respect to i5 = 7r/4. This 
is corroborated by exact analytical calculations based on the diagonalisation of Lj{x), which lead to a = 1/2 for 
A = 2, a == (55 - 12 sm{2S) - cos(4J))/36 for A^ = 4, and a = (303 - 110 sm{25) - 3 cos(4(5))/100 for A^ = 6 (in units 
7IG1 — G2p/(?'^). The plot shows that all numerical data can be very well fitted by a = A + Bsm{2S) + Gcos(45). 



From eq. (21) we know that A + C has to scale as A^^ for sufficiently large A 
which we have data, and C appears to be negligible 
rapidly than N'^ (a fit in the range A^ 



Supplementary Figure SI shows that 
, 18 gives a power law A^'^*) 



Both B and C are negative for all A^ for 
B increases even more 



But B has to cross over to a power law A^^ 
with p < 2, unless other Fourier components start contributing significantly. Otherwise, a would become negative for 
(5 > 0. This indicates that for large A^ the scaling of a is in fact A^^ for all 6. 

In summary, our method still works, even if the product state \iPq) is not prepared perfectly. One has the choice 
to start measurement immediately after state preparation, which gives an additional background of order A^, or to 
wait a time of the order of a few I/7 after preparation of the initial state, until no more photons leave the cavity 
through the mirror, with no additional background penalty. In both cases the scaling of the uncertainty of SL/L in 
a subsequent measurement of a small change of L is still ex 1/A^ . 
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Supplementary Figure SI: Effects of imperfect state preparation, (a) Scaling of collective photon escape rate a = 
— {Jz{x,0)) (in units 7IG1 — G2p/g^) in the mixture of DFS states reached by relaxation from state \ipo), eq. (|S6|I . Data 
for S = So {5 = ■k/2 — So) denoted by squares and full lines (diamonds and dashed lines), respectively; So = 0, 7r/9, 27r/9, Stt/Q 
and in/9 in black, red, green, blue, and brown. Exact analytical results for (5 = 0, eq. (21) shown with black crosses. Full and 
dashed lines are guides to the eye only, (b) Dependence of a on imperfection parameter S. N = 2, 4,6,8,10,12,14,16,18 in black 
circles, red squares, green diamonds, blue triangles up, brown triangles left, grey triangles down, violet triangles right, cyan 
pluses, orange Xs, respectively. The full lines for A^=2,4,6 are exact analytical results. The dashed lines for A'^ = 8, . . . , 18 are 
fits to A + Bsm{2S) + Ccos(4(5). The inset shows the scaling of the coefficients A (black circles), —B (red squares), and — C 
(green diamonds) as function of A''. 
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